Isoscaling of projectile-like fragments 
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The isotopic and isotonic distributions of the projectile fragmentation products have been sim- 
ulated by a modified statistical abrasion-ablation (SAA) model and the isoscaling behavior of 
projectile-like fragments has been discussed. The isoscaling parameters a and /3 for hot fragments 
before evaporation and cold fragments after evaporation have been extracted, respectively. It looks 
that the evaporation has strong effect on a. For cold fragments, a monotonic increase of a and \/3\ 
with the increasing of Z and -/V is observed. The relation between isoscaling parameter and the 
change of isospin content is demonstrated. In addition, the disappearance of the isospin effect of 
projectile fragmentation is also discussed in the viewpoint of isoscaling parameter. 

PACS numbers: 25.70.Mn, 42.10.Pa 



Recently, isoscaling behavior has been extensively ob- 
served 0,12, 0,3- The scaling law relates ratios of isotope 
yields measured in two different nuclear reactions, 1 and 
2, R 2l {N,Z) = Y 2 {N,Z)/Y 1 (N,Z). In multifragmenta- 
tion events, such ratios are shown to obey an exponential 
dependence on the neutron number N or proton number 
Z of the isotopes characterized by three parameters a, f3 
and C [jj]: 

R 21 (N,Z) = ^^- = Cexp(aN + (3Z). (1) 

In grand-canonical limit, a = A/j, n /T and f3 = A/i z /T 
where Afi n and Afi z are the differences between the neu- 
tron and proton chemical potentials for two reactions, re- 
spectively. C is an overall normalization constant. This 
behavior is attributed to the difference of isospin asym- 
metry between two reaction systems in the similar nu- 
clear temperature. It is potential to probe the isospin 
dependent nuclear equation of state [a, @ by the stud- 
ies of isoscaling 0, 0, 0, 0, ■ So far, the isoscaling 
behavior has been studied experimentally and theoreti- 
cally for different reaction mechanisms. However, most 
studies focus on the isoscaling behaviors for light parti- 
cles. A few studies on the heavy projectile-like residues 
in deep elastic collisions and fission fragments have been 
reported H E3 GH 

Il2l Il3j. In this paper, we concen- 
trate our attentions on the isoscaling features in projec- 
tile fragmentation and discuss the effect of evaporation 
on the isoscaling parameters a and (3 in a framework of 
statistical abrasion-ablation (SAA) model flij . 

Projectile fragmentation is one of the most important 
methods to produce beams of extreme neutron-rich or 
proton-rich nuclei and has been widely used to study nu- 
clear reactions induced by heavy ions at intermediate and 
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high energies. Various physical models for projectile frag- 
mentation have been developed and the reaction mecha- 
nisms of heavy-ion collisions have been investigated ex- 
tensively. For instance, an empirical parameterization, 
like EPAX/EPAX II formula presented by Summerer et 
al. |l5l Il6j | , can predict the fragment production cross- 
section well. However, due to its poor physical founda- 
tion, the empirical parameterization may not be appli- 
cable in extrapolation. In addition, statistical abrasion- 
ablation model can describe the isotopic distribution well 
[l4| . In the SAA model , the nuclear reaction is described 
in two stages which occur in two distinctly different time 
scales. The first abrasion stage is fragmentation reac- 
tion which describes the production of the prefragment 
with certain amount excitation energy through the in- 
dependent nucleon-nucleon collisions in the overlap zone 
of the colliding nuclei. The collisions are described by a 
picture of interacting tubes. Assuming a binomial distri- 
bution for the absorbed projectile neutrons and protons 
in the interaction of a specific pair of tubes, the distri- 
butions of the total abraded neutrons and protons are 
determined. For an infinitesimal tube in the projectile, 
the transmission probabilities for neutrons (protons) at 
a given impact parameter b are calculated by [T^| 

t k (r-b) = exp{-[D^(r-b)a nk +D^(r-b)a pk }}, (2) 

where D is thickness function of the target, which is 
normalized by J d 2 rl? n — N^- and J cPrDp = Z^ with 
T T 

N and Z referring to the neutron and proton number 
in the target respectively, the vectors r and b are defined 
in the plane perpendicular to beam, and cr k ' k is the free 
nucleon-nucleon cross sections (k', k=n for neutron and 
k', k=p for proton). The thickness function of the target 
is given by 

/-f-oc 
dzp k ((r 2 + z 2 ) 1/2 ), (3) 
-00 

with p k being the neutron (proton) density distribution 
of the target. So the average abraded mass at a given 
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FIG. 1: Yield ratios R,2i(N,Z) of pre-fragments (also called 
hot fragments in this work) from the reactions of 40 Ar + 9 Be 
and 3e Ar + 9 Be at 60 MeV/nucleon versus N (top panel) or 
Z (bottom panel). The lines represent the exponential fits. 



FIG. 2: Yield ratios R,2i(N,Z) of cold fragments from the 
reactions of 40 Ar + 9 Be and 36 Ar + 9 Be at 60 MeV/nucleon 
versus N (top panel) or Z (bottom panel). The lines represent 
the exponential fits. 



impact parameter b is calculated by the expression 

(AA(b))= jd?rDl(r)[l-tn(r-b)} 

+ Jcf r D^(r)[l-t p (r-b)}. { 1 

The excitation energy of projectile spectator is estimated 
by a simple relation of E* — 13.3(AA(b)) MeV where 
13.3 is a mean excitation energy due to an abraded 
nucleon from the initial projectile |l7^ . In the second 
evaporation stage the system reorganizes due to excita- 
tion, which means it deexcites and thermalizes by the 
cascade evaporation of light particles. By inducing in- 
medium nucleon-nucleon cross section and optimizing 
computational method proposed by our group given in 

Ref. EHUJEl, 

it can give a good agreement with 
the isotopic distribution at both high and intermediate 
energies involving neutron-rich or proton-rich nuclei over 
a wide energy range [l8l Hij. l2lj . The isospin effect and its 
disappearance in projectile fragmentation for 36 ' 40 Ar at 
intermediate energies have been predicted by this model 
and conformed by experimental data 19]. 

In order to study the isoscaling effect in projectile 
fragmentation, the reactions of 40 / 36 Ar + 9 Be at 60 
MeV/nucleon were simulated by SAA model. We extract 
the yield ratio R2i(N, Z) using the convention that index 
2 refers to the more neutron-rich system ( 40 Ar+ 9 Be) and 
index 1 to the less neutron-rich one ( 3e Ar+ 9 Be). Figure 
n shows the yield ratios i?2i(A, Z) of hot projectile-like 
fragments (PLFs) as a function of neutron number N 



for selected isotopes (upper panel) and proton number Z 
for selected isotones (lower panel) for the 40 / 36 Ar + 9 Be 
reactions in the semi-log plots. In figure |3 the corre- 
sponding ratios for the cold projectile-like fragments are 
shown in the semi-log plots. In the figures, the differ- 
ent isotopes and isotones are shown by alternating filled 
and open symbols: even isotopes and isotones are shown 
with filled symbols while odd ones are shown with open 
symbols. 

From the upper panels of Fig. (hot fragments) and 
Fig. (cold fragments), we observe that the ratio for 
each isotope Z (from 10 to 17) exhibits a remarkable ex- 
ponential behavior. For each isotope (Z), an exponential 
function form C exp(aiV) was used to fit the calculation 
points and the parameters a are shown in the top panel of 
Fig. [31for the selected Z. Analogous behavior is observed 
in the lower panels of Figs. ^and|21 for each isotone (A), 
an exponential function form C"exp(/?Z) was used to fit 
the calculation points and the parameters (3 are shown 
in the lower panel of Fig. [3] for the selected A. 

The model predicts that the fragment isotopic distri- 
bution at a fixed atomic number Z shifts towards the 
neutron-rich side for the neutron-rich projectile. In Fig. 
[3J we present the slope parameters a (upper panel) and 
\(3 1 (lower panel) of the exponential fits obtained as de- 
scribed for Figs, n and [21 as a function of Z and N, 
respectively. In the upper panel of Fig. [21 a of the pre- 
fragments for each element show a little bit decreasing 
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FIG. 3: Isoscaling parameters a as a function of Z (top panel) 
and /3 as a function of N (lower panel) for projectile fragmen- 
tation from the reactions 40 Ar + 9 Be and 36 Ar + 9 Be at 60 
MeV/nucleon. Solid symbols represent the pre- fragments and 
open symbols represent final fragments. 



trend with the increasing of Z. But, a monotonic increase 
of |/? with the increasing of N is observed. For cold frag- 
ments which are survived after evaporation (open sym- 
bols of Fig.|3), a monotonic increase of a and \/3\ with the 
increasing of Z and N, respectively, is observed. From 
the lower panel of Fig. [3J \(3\ always increase with the 
increasing of N for hot fragments or cold fragments and 
\[3 1 of hot fragments tends to drop when evaporation is 
taken into account. 

However, some cautions should be reminded. Since 
the first abrasion stage is essentially a fast geometrical 
abrasion stage, no equilibrium can be expected for the 
pre- fragments. These fragments are only intermediate 
stage products and will decay due to their excitation. In 
the strict sense, there is no isoscaling behavior for this 
pre-fragments due to the lack of equilibrium. However, 
a similar analysis as the isoscaling can be done since 
there exists different isotopic/isotonic distributions be- 
tween two systems. In this work, this kind of analysis is 
useful to investigate the evaporation effect on the isoscal- 
ing parameters of PLF fragments Q . 

In the evaporation stage, neutrons are emitted with the 
most important probability, the final cold fragments are 
more symmetric than their hot ancestors, which reflects 
a strong evaporation effect on projectile-like fragments 
p2l |23|. However, the charged particle emission is less 
important. Since a is the isoscaling parameter in fixed 



Z and P is the isoscaling parameter in fixed N, so the 
neutron evaporation has the strongest effect on a while 
charged particle evaporation has less effect on (3. That 
results in an obvious change of a parameter while a weak 
change of [3 due to the evaporation. 

In order to further discuss the reason of the dramatic 
change of the isoscaling parameters, we plot (N)/Z ver- 
sus Z in the upper panel of Fig. 21 for the fragments be- 
fore and after evaporation. As known, (N) jZ represents 
the isospin component. From the figure, the hot PLFs 
which are close to initial projectile basically remain larger 
isospin content as the initial projectile while lighter hot 
PLFs have much larger (N)/Z due to stronger dissipation 
and nucleon exchange. However, this high isospin content 
can not keep alive since the prefragments are excited as 
illustrated the above section. They will cool themselves 
by the light particle evaporation, mostly by the neutron 
emission. After the evaporation process, (N)/Z tends to 
1.1 in lower Z and bifurcates in higher Z for Ar and 
36 Ar systems. More interestingly, we note that the dif- 
ferences of (N)/Z between two systems for hot and cold 
fragments as a function of Z as shown in the lower panel 
of Fig. Consistently, this behavior is similar to a as a 
function of Z for hot and cold PLFs as shown in Fig. 
In this context, we could say the difference of isospin con- 
tent (i.e. the neutron-to-proton ratio: (N) jZ) is a more 
simple measurement of isoscaling parameter a. Similar 
behavior has been found for 

In the past years, the isospin dependence of various 
physical quantities has been reported [jj 0, 0, |25|, |2(| 
|27| . The positive slopes in the upper panels of Figs. 
and |3 indicate that neutron-rich fragments are more ef- 
ficiently produced, as expected, from the more neutron- 
rich systems. In other words, the fragment isotopic dis- 
tribution at a fixed atomic number Z shifts towards the 
neutron rich side for the neutron-rich projectile. This 
isospin effect will decrease with decreasing the fragment 
atomic number Z and disappear when (Z pro j — Z) /Z pro j 
becomes larger ,19j. The quantity of {Z pro j — Z)/Z pro j 
can be related to excitation energy or temperature of 
the deformed projectile after the abrasion stage. When 
this quantity becomes smaller, the deformed projectile 
becomes hotter or highly excited. In this context, we 
can say that the isospin effect tends to fade out with the 
increasing of the excitation energy of the projectile-like 
system. 

In summary, the isoscaling of projectile-like fragments 
from 40 / 3( >Ar + 9 Be at 60 MeV/nucleon has been studied 
by a modified statistical abrasion-ablation model. SAA 
model can simulate well the isotopic distribution of PLFs 
of the systems. The isoscaling parameters a and (3 are 
extracted for the prefragments and final fragments, a 
has a dramatic change after evaporation. For hot PLFs, 
a shows a little bit stronger isospin effect for those highly 
excited small PLFs. However, this kind of stronger 
isospin effect for smaller PLFs can not survive due to 
much more neutron evaporation from these highly ex- 
cited pre-PLFs. Finally a shows weaker for these smaller 
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FIG. 4: (a) (N)/Z for hot and cold fragments for 40 Ar and 
36 Ar systems. The dotted line is initial N/Z for 40 Ar; (b) 
The differences of (N)/Z between 40 Ar and 3e Ar for the hot 
and cold fragments systems. The symbols are illustrated in 
insert. 



PLFs while stronger for PLFs close to initial projectile. 
A monotonic increase of a ( |/3| ) with increase Z ( N 
) is observed. This behavior can be easily explained by 
the difference of the isospin content of hot or cold PLFs 
between two systems as shown in Fig. 0] In this sense, 
isospin content can take a simple "isometre" as isoscaling 
parameter. In addition, the disappearance of the isospin 
effect of projectile fragmentation is also discussed in the 
viewpoint of isoscaling parameter. 
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